Breakup of a Dimer: A New Approach to Localization Transition 

Ignacio Gomezld and Indubala I. Satijaa 

Department of Physics, George Mason University, Fairfax, VA 22030 

(February 1, 2008) 



00 



> 

o 

o 

0\ 
-(— > 



X3 
O 

o 



X 



Within the framework of tight binding models, aperiodic 
systems are mapped to a renormahzed lattice with a dimer 
defect. In models exhibiting metal-insulator transition, the 
dimer acts like a resonant cavity and explains the existence 
of the ballistic transport in the system. The localization in 
the model can be attributed to the vanishing of the coupling 
between the two sites of the dimer. Our approach unifies 
Anderson transition and resonance transition and provides a 
new formulation to understand localization and its absence in 
aperiodic systems. 

PACS numbers: 72.15.Rn+72.15-v 



The existence of metal-insulator transitions due to 
quantum interference in correlated disordered systems is 
a fascinating subject. The origin of this resonant tran- 
sition can be understood from the textbook examples of 
the quantum mechanics of a barrier in the continuum 
problem and from a simple dimer model |^ in the case 
of a discrete lattice. The resonance transition where the 
metallic phase with ballistic transport is due to zero re- 
flectance at the two sites of the dimer provides a simple 
mechanism to understand localization and its absence in 
systems with short range correlation. Recently, this tran- 
sition was verified in experiments in superlattices Q . 

In this letter, we present a universal approach to un- 
derstand metal-insulator transitions by unifying Ander- 
son transition with the resonance transition. Our for- 
mulation is applicable for any aperiodic system with re- 
flection symmetry described by a nearest-neighbor tight 
binding model (TBM). In this paper we will confine our- 
selves to sinusoidal potentials which are further modu- 
lated by a Gaussian profile. The purpose of the Gaussian 
modulation is two-fold: as explained below, it provides 
a useful means to motivate and illustrate our ideas even 
though our results and conclusions are valid in the limit 
where the width of the Gaussian goes to oo. Secondly, it 
facilitates the study of the pure Gaussian systems that 
have been the subject of recent studies [Q due to its pos- 
sible application as efficient energy band-pass filters in 
semiconducting superlattices. Finally, the case of sinu- 
soidal potential further modulated by a Gaussian profile 
includes the famous Harper equation exhibiting Ander- 
son localization as a limiting case. Q 

The model system under consideration here is a TBM 
describing an eigenvalue problem with energy E, 



The Em = E ~ e„i is the diagonal term containing the 
aperiodic onsite energy e„i. This model also describes the 
Schrodinger equation for an array of (5-function Kronig- 
Penny potential barriers, 

-^ip{x)+'^eraS{x-ma)'ilj{x)=Eil;{x). (2) 

This is due to the fact that the Poincare map m associ- 
ated with the model is a TBM with _E,„ = e„i ''"^ ' — 
cos{K). Here K is the Bloch vector related to the energy 
as ^ = h^K^/2M. 

We choose em to be sinusoidal potential which is fur- 
ther modulated by a Gaussian, 
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Here A is the strength of the potential, a is the width of 
the Gaussian profile and /i is an irrational number chosen 
to be the inverse golden-mean ^ — ^^^ ■ In this paper, 
we will describe our results for the pure Gaussian model 
(/i = 0), the Harper equation {a — > cxd) and also for the 6- 
function barriers Eq. (g) The last case closely resembles 
the recent study of Gaussian modulated Kronig-Penny 
model g and is also currently under experimental inves- 
tigation in superlattices. g]. 

The basic idea underlying our approach can be under- 
stood for any localized defect with finite spatial extent. 
For aperiodic systems such as Harper equation where the 
aperiodicity exists throughout the lattice, we introduce 
Gaussian modulation so that the system can be viewed 
as a lattice with a localized defect. However, as shown 
below, the ct — > cx) limit is well defined and therefore 
Gaussian modulation , although not needed, is useful in 
understanding the decimation scheme described below. 

We envision perfectly transmitting phase in all aperi- 
odic systems to be described by Bloch states on some 
" renormalized" lattice. In aperiodic systems where the 
defects are spatially confined to only some parts of the 
lattice, Bloch wave amplitudes are attenuated or ampli- 
fied due to various quantum interferences only in the 
neighborhood of these defects. ( See figure 1) We elimi- 
nate such sites from the model by decimating them. The 
resulting renormalized model will have solutions that are 
Bloch waves at all sites in the metallic phase. We would 
like to emphasize that we decimate all sites which is in 
contrast to previous use of decimation for random sys- 
tems where every other site is decimated. However, the 
basic idea of decimating all sites is similar in spirit to that 



of Fibonacci decimation |Q] used in quasiperiodic systems 
where one eliminates all but Fibonacci sites so that the 
renormalized lattice may exhibit translational invariance 
in Fibonacci space. 

Figure 2 outhnes the decimation process. We begin by 
eliminating the central site mo of the symmetric defect. 
This leads to a renormalized model consisting of a dimer 
with coupling between the two sites as 71 = I/Eq and the 
onsite energy Ei = Ei — ji. We now begin the iterative 
process by decimating this dimer. At the n*'' step, we 
obtain a new lattice with renormalized dimer with onsite 
energies denoted as En+i and the renormalized coupling 
between the 2-sites of the dimer as 7n+i. This iterative 
decimation scheme results in a two-dimensional driven 
map for j — E where the En containing diagonal disorder 
of the bare model provides the driving term, 
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It should be noted that the parameters of the model 
are included in En+i- This 2-dimensional map contains 
the complexity of various interference effects within the 
Gaussian profile, manifesting itself in the energy depen- 
dent coupling and onsite potential. We would like to 
emphasize again that the map has a well defined limit 
for a equal to 00. Therefore, Gaussian profile is not nec- 
essary in obtaining the mapping of the aperiodic model 
to the dimer model. The dynamics of the map does de- 
pend upon a: the finite value of a provides damping in 
the map with the consequence that the renormalization 
group (RG) flow settles on attractors while, in the limit 
of cr — > 00, there are no attractors in the map. 

The usefulness of this map emerges from the fact that 
the changes in the parameters A are reflected in the sig- 
nificant changes in the trajectories of the map for both 
finite as well as for the infinite value of a. Figure 3 shows 
how the metal-insulator transition in the pure Gaussian 
case manifests itself in terms of the variation in the RG 
attractor as A is varied. In the subcritical phase, where 
the lattice has a transmission coefficient T equal to unity, 
a symmetric period-2 limit cycle describes the RG flow 
for the renormalized coupling 7 while the dynamics of E 
is governed by a fixed point. Away from the transition 
point, the RG attractor exhibits very regular oscillatory 
pattern consisting of symmetric lobes and divergences 
which is periodic in A with periodicity equal to -. The 
fact that the total number of lobes or divergences is equal 
to a suggests that each lobe owes its existence to a partic- 
ular site within the Gaussian profile where the sites near 
the center contributing at smaller values of A. However, 
the quantitative understanding of almost equally spaced 
lobes and its physical significance remains elluded to us 
at present. 

As we approach the transition point, the symmetric 
2-cycle of 7 looses its symmetry, degenerates to a fixed 



point (seen in the figure 3 where the lobes cross) and then 
continues as an asymmetrical period-2 attractor. On the 
other hand, the fixed point describing the E values be- 
comes a 2-cycle. At the transition, 7 approaches zero 
with a power-law decay resulting in a broken dimer. The 
localized phase is characterized by an exponential vanish- 
ing of the the coupling with a characteristic length which 
is found to be related to the localization length ^, of the 
localized wave function, 7„ « exp — {2n/^). It turns out 
that in the localized phase, E continues to be described 
by a period-2 fixed point exhibiting divergences where 
the spacing between two successive divergences increases 
as A increases. 

The localization transition discussed above is a reso- 
nant transition where the metallic phase is described by 
Bloch wave solutions which undergo real phase shifts as 
they encounter the dimer defect. In the localized phase, 
these phase shifts become imaginary. Imposing this con- 
dition on the solutions of the renormalized system de- 
termines the condition for the perfect transmission. We 
express this condition in terms of a function / defined as 
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Figure 4 shows the variation in /„ and the transmission 
coefficient T for the Kronig-Penny model Eq. (||) as the 
energy E is varied. We see that the condition for the 
vanishing of this function coincides with the condition 
for perfect transmission. We would like to point out that 
the transmission coefficient was calculated in a rather 
simple way by using the renormalized dimer [n2|. This 
requires multiplying only two transfer matrices in con- 
trast to usual calculations where one needs to multiply 
all transfer matrices on an aperiodic lattice. 

The existence of a band of conducting states in Gaus- 
sian modulated lattices as discussed above ( which is con- 
sistent with the previous related study g] ) is a very desir- 
able feature of a lattice making it a useful filter. Attempts 
are underway |g| to make such filters using superlattices. 
This should be contrasted with the dimer-type defects 
where the metal-insulator transition is obtained by fine 
tuning the parameter to obtain the resonant energy close 
to the Fermi energy. We would like to point out that the 
metal-insulator transition and a band of perfectly trans- 
mitting states in Gaussian system can also be understood 
by heuristic arguments of asymptotic property of "con- 
stantcy" of the potential. Due to the finite width of the 
potential, the model has a constant potential asymptot- 
ically. This property is crucial in determining the E, 
which is the global property lO]. This argument leads 
to extended states if ii^ > 2 — 2A and localized states for 
E < 2 — 2A. This condition is found to be true in our 
numerical simulations for large values of a. 

We next discuss the case of the pure Harper equation 
obtained in the limit cr ^ 00 in the driving term for the 
two-dimensional map. The sub-critical phase is not an 



attractor. As A increases, the RG trajectories become 
more and more complex ( see figure 5) [ [L3[ and eventu- 
ally collapse to a vertical line corresponding to 7 -^ at 
the onset to localization transition. The localized phase 
is again characterized by an exponentially decaying cou- 
pling of the dimer with length scale which is equal to 
half the localization of the Harper equation. The metal- 
lic phase in Harper equation is thus described in terms of 
resonace due to the dimer and the Anderson localization 
is due to the breaking of this dimer. 

Recently, the existence of extended states in super- 
critical Frenkel-Konterova(FK) model and in Fibonacci 
lattices was shown to be due to dimer-type correlations 
|10| using decimation schemes that were model depen- 
dent. The novel aspect of the work described here is 
the universal nature of our approach: the 2-dimcnsional 
map Eq. (Q) can be used to study localization or its ab- 
sence in systems with short range correlations such as 
dimer defects, FK and Fibonacci models with long range 
correlations as well as in Harper equation. Our most 
important result is the picture of the localized phase as 
a phase with the broken dimer. Although the present 
framework is developed only for the TBM type systems, 
it is possible to extend our method to study localiza- 
tion in 2-dimensional models as well as for systems with 
long range interaction. The new methodology developed 
here will have applications in many other areas includ- 
ing dynamical localization in kicked rotors WM as well as 
the transition to strange nonchaotic attractors(SNAs) in 
quasiperiodically driven maps [ p^ , p^ . 
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FIG. 1. The wave function for sinusoidally modulated 
Gaussian case. We have Bloch wave with K = 1/2 outside 
the Gaussian profile. The points in the wave function are not 
joined for clarity. Here A = 0.5 and a = 1000. 



FIG. 2. The Gaussian modulated lattice and the iterative 
decimation process. 



FIG. 3. The variation in the asymptotic values of the 
renormalized coupling 7„ and En for a = 100. There are 
few isolated points in the localized phase near the transition 
where 7„ may diverge. Note that tanh is used to display the 
behavior at ±00. 



FIG. 4. The band of resonant energies corresponding to 
perfect transmission for (5-function potential Kronig-Penny 
model. Here A = 10 and a = 100. 
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FIG. 5. RG trajectory for A ■- 
Harper equation. 
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